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Molecular  beams  of  hydroqen  bonded  water  dimer,  generated  in  a 
supersonic  nozzle,  have  been  studied  using  electric  resonance 
spectroscopy.  Radio,  frequency,  and  microwavef^  transitions  have 
been  observed  in  (H2'^%)2.  (02^60)2,  and  (H2^®0)2.  Transitions 

arising  from  both  pure  rotation  and  rotation-tunneling  occur. 
The  pure  rotational  transitions  have  been  fit  to  a rigid  rotor 
model  to  obtain  structural  information.  Information  on  the 
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ABSTRACT 

Molecular  beams  of  hydrogen  bonded  water  dimer,  generated  in 

a supersonic  nozzle,  have  been  studied  using  electric  resonance 

spectroscopy.  Radio  frequency  and  microwave  transitions  have  been 
X6  16  18 

observed  in  {H2  0)2'  ^^2  2 ' ^^2  ^^2’  Transitions  arising 

from  both  pure  rotation  and  rotation-tunneling  occur.  The  pure  ro- 
tational transitions  have  been  fit  to  a rigid  rotor  model  to  obtain 
structural  information.  Information  on  the  relative  orientation  of 
the  two  monomer  units  is  also  contained  in  the  electric  dipole 
moment  component  along  the  A inertial  axis,  y , which  is  obtained 

Si 

from  Stark  effect  measurements.  The  resultant  structure  is  that  of  a 
" trans-linear"  complex  with  an  oxygen-oxygen  distance,  Rqq^  of  2.98(1) 

O 

A,  the  proton  accepting  water  axis  is  58(6)°  with  respect  to  Rqqc  and 
the  proton  donating  water  axis  at  -51(6)°  with  respect  to  R^^.  This 
structure  is  consistent  with  a linear  hydrogen  bond  and  the  proton 
acceptor  tetrahedrally  oriented  to  the  hydrogen  bond.  The  limits  of 
uncertainty  are  wholly  model  dependent  and  are  believed  to  cover 
variations  from  the  zero-point  vibrational  structure  observed  to  the 
equilibrium  structure.  y^  shows  strong  dependence  on  J and  K and  is 
about  2.6D.  Centrifugal  distortion  constants  have  been  interpreted  in 
terms  of  the  monomer-monomer  stretching  frequency  and  give  w = 150  cm 
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INTRODUCTION 

The  structure  of  liquids  in  general,  and  the  structure  of  liquid 
water  in  particular,  have  been  the  subject  of  a vast  amount  of  re- 
search. However,  the  complexity  of  the  problem  is  so  great  that  only 
initial  successes  have  been  achieved  and  liquid  structure  will  remain 
an  active  research  field.  The  traditional  approach  to  the  structure 
of  liquid  water  is  to  consider  a macroscopic  sample.  Experimental 
studies  have  applied  many  spectroscopic  and  diffraction  techniques  to 
water,  as  well  as  most  classical  physical  and  chemical  measurement 
methods.^  A great  variety  of  theoretical  models  have  been  applied  to 

liquid  water  ranging  from  mixtures  of  ice-like  crystallites  to  molecular 

2 

dynamics  calculations. 

In  the  past  several  years  there  has  been  a growing  emphasis  on 
approaching  the  liquid  water  problem  from  a microscopic  point  of  view. 
The  focus  of  this  approach  is  the  hydrogen  bond  between  two,  or  a small 
number,  of  isolated  water  molecules.  Thus,  the  species  studied  are  hy- 
drogen bonded  water  dimers,  trimers,  etc.  A detailed  knowledge  of 
these  isola::ed  systems  will  provide  a more  complete  understanding  of  the 
hydrogen  bond  and  of  intermolecular  potential  functions  that  dominate 
the  structure  of  liquid  water.  Such  studies  of  water  and  other  small 
hydrogen  bonded  systems^  will  lead  to  a much  more  general  knolwedge  of 
all  hydrogen  bonding.  The  detailed  study  of  systems  such  as  water  dimer 
should,  therefore,  contribute  to  a wide  range  of  chemical  problems  in- 
cluding structure  of  condensed  phases,  solvation,  intramolecular  s>truc- 
ture,  and  many  areas  of  biochemistry  involving  hydrogen  bonding. 
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A very  large  number  of  ^ initio  and  semi-empirical  quantum  cal- 
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culations  have  been  performed  on  water  dimer.  As  these  calculations 
have  become  more  sophisticated  and  contain  fewer  restrictions,  a 
consensus  on  the  theoretical  structure  of  water  dimer  has  arisen.  The 
lowest  energy  equilibrium  structure  calculated  is  the  " trans-linear" 
form.  This  structure  contains  a plane  of  symmetry  and  a linear,  or 
near  linear,  hydrogen  bond.  A much  smaller  number  of  calculations 
have  been  performed  on  the  trimer  and  larger  polymers  and  no  structural 
consensus  exists.^  In  particular,  there  is  no  agreement  as  to  whether 
the  trimer  is  an  open  structure  containing  two  hydrogen  bonds,  or  a 
cyclic  entity  with  three  hydrogen  bonds. 

From  an  experimental  point  of  view,  water  dimers  were  first  postu- 
lated to  explain  PVT  observations  on  water  vapor. ^ Since  the  initial 
concept  of  isolated  hydrogen  bonded  species,  there  have  been  many  direct 

and  indirect  experimental  investigations  of  these  molecules.  Water 
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dimers  have  been  directly  observed  using  matrix  isolation  infrared  ' 
and  mass  spectrometer  techniques . ^ The  initial  infrared  work  of 
Pimentel  and  coworkers^  was  interpreted  in  terms  of  a symmetric  structure 
containing  two  strongly  bent  hydrogen  bonds.  This  conclusion  was  based 

p 

on  the  observed  number  of  normal  modes.  Later  infrared  work  has  ob- 
tained more  data  and  concluded  the  dimer  structure  was  more  open  with  a 
single  hydrogen  bond,  in  agreement  with  theoretical  predictions.  While 
mass  spectroscopy  is  normally  insensitive  to  structural  details,  an 
electric  deflection  study  of  water  dimer  showed  the  symmetric  structure 
of  Pimentel  js  not  correct. Preliminary  communcations  on  the  microwave 


spectroscopic;  work  reported  here  showed 


water  dimer  to  have  the 
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"trans-linear"  structure.  Much  less  experimental  information  on 
the  structure  of  water  trimer  is  available.  The  electric  deflection 
observation^^  of  a very  small  permanent  electric  dipole  for  the  trimer 
is,  however,  a strong  piece  of  evidence  favoring  the  three  hydrogen 
bond  cyclic  structure. 

This  paper  presents  the  details  of  a molecular  beam  electric 
resonance  study  of  the  radio  frequency  and  microwave  spectra  of  water 
dimer.  The  spectroscopic  observations  are  interpreted  in  terms  of 
a quantita .live  geometric  structure  for  the  ground  vibrational  state 
of  water  dimer.  Electric  dipole  moment  measurements  are  also  presented 
The  nonrigid  nature  of  water  dimer  adds  several  complexities  to  its 
rotational  spectra.  The  preceding  paper^^^  discusses  the  theoretical 
aspects  of  this  nonrigid  behavior. 

EXPERIMENTAL 

Apparatus 

A molecular  beam  of  hydrogen  bonded  water  dimer  was  formed  by  ex- 
panding water  vapor  through  a simple  pinhole  nozzle. Nozzles  of 
different  diameters  and  thicknesses  were  investigated  and  optimum  dimer 
beams  were  produced  using  0.1  mm  diameter  and  0.13  mm  thickness.  The 
nozzle  was  fabricated  by  using  a jeweler's  drill  to  place  the  hole  in 
nickel  foil  of  chosen  thickness. The  foil  was  sealed  to  a conven- 
tional beam  source  using  an  indium  gasket.  A skimmer  was  not  used.  How 
ever,  a liquid  nitrogen  cooled  cryopump  with  a one  centimeter  aperature 
for  the  beam  was  located  about  two  centimeters  in  front  of  the  nozzle. 
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As  has  been  reported  previously,  ' this  type  of  beam  source 
can  produce  large  water  clusters  if  relatively  high  pressures  of  water 
vapor  are  used.  At  one  atmosphere  source  pressure,  all  (H^O)^  clusters 
within  the  m/e  = 350  mass  limit  of  the  molecular  beam  mass  spectrometer 
detector  have  been  observed. Optimum  production  of  dimers  was  ob- 
served at  roughly  100  torr  source  pressure  and  80 °C  source  temperature. 
The  source  pressure  has  not  been  directly  measured  but  was  adjusted 
to  maximize  the  dimer  concentration  in  the  beam.  The  dimer  content  of 
the  beam  was  monitored  at  the  ion  peak  of  the  mass  spectrometer 

beam  detector.  Since  higher  polymers  also  contribute  fragment 

ions,  dimer  concentration  was  judged  by  observing  only  those  molecules 
which  were  focused  around  a beam  stop  by  the  state  selecting  fields  of 
the  electric  resonance  spectrometer.  Since  the  higher  polymers  are 
essentiallj'  nonpolar^^  they  were  not  focused  into  the  detector.  Under 
standard  operating  conditions  the  ion  peak  from  water  dimer  was  ap- 

proximately 10%  as  intense  as  the  H20^  peak  from  water  monomer. 

Because  it  was  necessary  to  use  rather  low  source  pressure  to  avoid 
higher  polymer  formation,  the  absolute  intensity  of  the  dimer  beam  was 
small.  The  low  source  pressure  also  meant  the  nozzle  was  not  operating 
tar  into  the  supersonic  region  and  little  cooling  of  rotational  degrees 
of  freedom  v/as  observed. Therefore,  a relatively  weak  molecular  beam 
with  a large  partition  function  was  used  and  individual  rotational 
transitions  were  observed  with  relatively  poor  signal  to  noise  ratio. 

2 beams  were  formed  from  normal  water^^  and  were  monitoi’od  at 
m/e  * 19.  were  generated  from  99%  enriched  D^O  vapor  and 

were  detected  at  D^O  , m/e  =22.  A 93%  °0  enriched  sample  of  H^O, 


-6- 


kindly  provided  by  Dr.  Pinchas  of  Weizman  Institute,  was  used  to  pro- 

18  18  + 
duce  (H2  0)2  species.  These  molecules  were  observed  at  O at 

m/e  = 21.  One  day  of  operation  required  several  grams  of  water  so 

fairly  large  amounts  of  material  were  required.  Most  of  the  water  is 

1 8 

trapped  on  the  cryopump  and  the  small  H2  O sample  was  recycled. 

The  molecular  beam  electric  resonance  spectrometer  used  in  this 

1 7 

study  has  been  previously  described.  Modifications  not  discussed 
earlier  include  the  molecular  beam  source,  microwave  electronics,  and 
signal  averaging  capability.  Microwave  radiation,  generated  by  reflex 
klystron  oscillators,  was  introduced  into  the  resonance  region  of  the 
spectrometer  by  a waveguide  which  extends  up  to  the  electrodes  which 
provide  the  uniform  Stark  field.  These  electrodes  (the  C-field)  allow 
the  microwave  radiation  to  propagate  from  the  waveguide  to  the  beam 
molecules.  In  a nearly  prolate  rotor  with  the  K quantum  number  esqual 
to  zero  it  is  desirable  to  observe  AM=0  transitions  and  in  this  case 
the  broad  wall  of  the  waveguide  is  parallel  to  the  beam  axis  to  make 
the  radiation  and  static  fields  parallel. 

It  is  not  obvious  that  rotational  transitions  in  symmetric  top 
molecules  with  K>0  can  be  observed  in  electric  resonance  spectroscopy. 
The  requirement  that  the  Stark  coefficient  of  the  initial  and  final 
level  be  greater  and  less  than  zero,  respectively,  certainly  cannot  be 
met  with  parallel  radiation.  Water  dimer  a-type  rotational  transitions 
with  K^2  are  essentially  identical  to  those  of  a perfect  symmetric 
top  since  the  Stark  shift  is  typically  larger  than  the  asymmetry  split- 
ting. Such  transitions  have  been  observed  with  perpendicular  micro- 
wave  radiation.  The  feasibility  of  this  type  of  transition  was  ascer- 


tained  by  observing  virtually  identical  rotational  transitions  in 

2 18 
the  02  0 vibrational  state  of  carbonyl  sulfide.  The  observation 

of  these  transitions  can  be  explained  in  two  ways.  It  is  possible 
to  observe  a transition  from  a level  with  a positive  Stark  coeffi- 
cient to  a level  with  essentially  zero  Stark  coefficient  and  thus 
the  transition  can  be  M = 1 ->•  M'  = 0 . It  is  also  possible  to  start 
with  M>0  and  achieve  M'<0 levels  by  inducing  multiple  AM  = 1 transi' 
tions,  each  transition  reducing  M by  one,  to  cause  a cascade  of 
AM>>1.  Since  observations  were  always  made  with  excess  microwave 
power,  the  latter  explanation  is  probably  important. 

Several  different  procedures  were  used  in  measuring  microwave 

transitions.  For  initial  searches  it  was  desirable  to  have  linewidths 

as  large  as  possible  and  this  was  achieved  by  frequency  modulating 

19 

the  radiation  source  with  random  noise.  Since  the  radiation  was  not 
detected,  this  technique  does  not  lower  the  sensitivity  of  the 
spectrometer  as  long  as  sufficient  power  per  unit  bandwidth  is  available. 
With  approximately  100  mw  of  microwave  power,  linewidths  in  excess  of 
1C  MHz  were  readily  achieved.  Convenient  modulation  is  obtained  by 
mechanically  switching  the  noise  voltage,  which  is  applied  to  the  kly- 
stron repeller,  on  and  off  at  a 10  Hz  rate.  Transition  frequencies  can 
be  measured  to  10  MHz  accuracy  using  a cavity  frequency  meter  and  1 MHz 
accuracy  war,  obtained  by  observing  zero  beating  with  a known  frequency 
during  the  half  modulation  cycle  when  the  noise  was  off.  For  higher 
resolution  measurements  it  was  necessary  to  phase  lock  the  klystron  to 
a stable  frec(uency  derived  from  a crystal  oscillator  and  frequency  syn- 
thesizer. The  linewidth  observed  is  about  30  KHz  FWHM  which  results  from 
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unresolved  hyperfine  structure  and  non-optimum  radiation  distribution 
in  the  resonance  region.  Initial  measurements  were  made  using  fre- 
quency modulation  accomplished  by  modulating  the  30  MHz  reference  oscil- 
lator of  the  phase  locked  loop. 

It  was  difficult  to  maintain  phase  lock  and  adequate  frequency 
deviation  and  still  sweep  the  microwave  frequency.  An  amplitude 
modulator  capable  of  100%  modulation  at  10  Hz  over  the  8-60  GHz  fre- 
quency region  is  necessary  to  simplify  measurements  of  this  kind.  A 
simple  mechanical  modulator  was  built  for  this  purpose.  A short  length 
of  X-band  waveguide  was  fitted  with  an  oscillating  vane  connected  to  a 
small  shaft  passing  through  the  center  of  the  narrow  wall  of  the  wave- 
guide. KTien  this  vane  is  parallel  to  the  broad  wall  in  the  middle  of 
the  waveguide  it  is  a perpendicular  to  the  microwave  E field  and  causes 
minimum  attenuation.  Rotation  of  the  vane  shaft  approximately  20  de- 
grees brings  the  edges  of  the  vane  in  contact  with  the  broad  walls  of 
the  guide,  blocking  microwave  transmission.  The  shaft  is  made  to 

oscillate  at  10  Hz  by  a galvanometer  drive  manufactured  for  optical 
20 

scanning  mirrors.  Cross  sectional  sketches  of  the  modulator  are  shown 
in  figure  one.  Using  tapered  transition  sections  this  modulator  can 
be  used  for  frequencies  throughout  the  usual  microwave  spectroscopic 
region.  With  this  modulator  phase  locked  klystrons  were  stable  and 
easily  swept.  Signal  averaging  was  accomplished  by  interfacing  a multi- 
channel analyzer  to  the  frequency  synthesizer  in  the  phase  locked  loop. 

Radio  frequency,  RF,  transitions  were  observed  by  generating 
suitable  radiation  fields  both  parallel  and  perpendicular  to  the  static 
field.  The  RF  source  is  a Hewlett-Packard  3330A  frequency  synthesizer. 
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The  frequency  range  is  extended  by  harmonic  generation  and  filtering 
using  a Hewlett-Packard  10511A  spectrum  generator  and  230B  amplifier. 

Radio  frequency  transitions  are  10-30  KHz  FWHM  due  to  unresolved  hyper- 
fine  structure.  Instrumental  linewidth  is  about  3 KHz. 

The  frequencies  reported  are  the  average  of  data  taken  sweeping 
frequency  up  and  down.  Most  measurements  were  made  using  a 30  second 
detector  time  constant  with  signal  to  noise  ratios  varying  between  5 
and  20.  Signal  averaging  was  used  to  improve  signal  to  noise  ratios 
for  some  lines,  particularly  for  Stark  effect  measurements.  Figure  2 
shows  signal  averaged  results  for  two  radio  frequency  transitions. 

The  assignment  of  many  microwave  lines  was  confirmed  by  RF-micro- 
wave  double  resonance  experiments.  These  experiments  were  most  fre- 
quently carried  out  by  fixing  the  modulated  RF  frequency  at  an  assigned 
transition  and  monitoring  the  intensity  of  this  RF  transition  as  the 
unmodulated  microwave  frequency  is  swept.  The  intensity  of  the  RF 
transition  is  observed  to  either  increase  or  decrease,  depending  on  the 
relative  energy  level  scheme,  when  both  transitions  occur  simultaneously 
and  involve  a comjnon  energy  level. 

Preliminary  Discussion 

All  of  the  spectroscopic  transitions  we  have  observed  in  water 
dimer  fall  into  one  of  two  groups.  Many  transitions  clearly  arise  from 
a nearly  prolate  symmetric  top  rotor  with  permanent  electric  dipole 
moment  component  parallel  to  the  A axis,  which  is  the  approximate  pro- 
late axis.  The  remaining  transitions  show  no  apparent  pattern  and 
will  not  fit  a rigid  rotor  model.  As  was  discussed  in  detail  in  the 

1 
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previous  paper, these  observations  are  consistent  with  hydrogen 
bonded  water  dimer  molecules  capable  of  undergoing  proton  tunneling 
motions  which  interchange  identical  nuclei.  The  majority  of  micro- 
wave  transitions  will  involve  both  rotational  and  tunneling  energies. 

It  was  shown  that  there  are  at  least  six  tunneling  splittings  of  un- 
known magnitude  and  the  resultant  spectrum  is  indeed  very  complicated. 
There  are,  however,  allowed  transitions  between  two  of  the  tunneling 
sublevels  which  only  involve  rotational  energies,  i.e.,  these  transi- 
tions connect  the  same  tunneling  sublevels  in  two  different  rotational 
states.  Tie  two  tunneling  sublevels  in  question  are  doubly  degenerate 
with  symmetry  and  E . Figures  3 and  4 of  the  previous  paper  show 
the  allowed  transitions  for  water  dimer. 

The  picture  of  water  dimer  undergoing  large  amplitude  tunneling 
motions  was  not  available  when  this  study  was  undertaken.  Rather, 
this  model  was  developed  to  explain  the  observed  spectroscopic  data. 
This  approach  was,  however,  greatly  facilitated  by  the  earlier  study 
of  hydrogen  fluoride  dimer  which  was  observed  to  exhibit  tunneling.^ 

The  frequencies  obtained  from  the  pure  rotational  transitions 

between  the  E states  can  be  fit  with  conventional  rigid  rotor  theory 

22 

using  effective  rotational  constants.  Water  dimer  is  necessarily  a 

slightly  as/mmetric  prolate  top  since  it  contains  only  two  heavy  atoms. 

The  A inertial  axis  will  be  nearly  coincident  with  the  oxygen-oxygen 

axis  for  any  structure  of  the  molecule.  The  A,  R,  and  C rotational 

constants  can  bt  estimated  from  approximate  structures  based  on  oxygen- 

23  24 

oxygen  distances  obtained  from  studies  of  ice  or  liquid  water,  or 

4 

structures  from  a^  xnitio  or  semi-empirical  calculations  can  be  used. 


i 
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5 3 3 

Rough  estimates  are  A*2xl0  , B=6xl0  , C*6xl0  MHz.  The  degree  of 
asymmetry  is  highly  dependent  on  B-C  which  is  determined  by  fine 
structural  details.  B-C  can  be  expected  to  be  in  the  range  of  0 to 
100  MHz.  Thus  the  molecule  is  strongly  prolate  and  only  slightly 
asymmetric.  The  vector  sum  of  the  monomer  dipole  moments  will 
guarantee  a large  component  of  electric  dipole  moment  parallel  to  the 
A axis.  The  perpendicular  component  of  the  moment  depends  critically 
on  the  relative  orientation  of  the  monomer  units  and  can  vary  from 

zero  to  greater  than  cwo  Debyes.  j 

With  tnese  estimates,  rigid  rotor  spectra  can  be  calculated.  In  j 

2 ^ 

the  synunetric  top  limit,  the  energy  levels  are  given  by  BJ ( J+1 ) + (A-B)  K j 

where  J and  K specify  the  total  angular  momentum  and  the  projection  | 

of  that  momentum  on  the  A axis.  Transitions  associated  with  u , the  j 

a 

dipole  mom»!nt  parallel  to  the  A axis,  are  referred  to  as  A type  transi- 
tions and  change  J by  one  unit  and  leave  K unchanged.  These  trar.si-  | 

1 

tions  occur  at  2B,4B,6B,  etc.  Perpendicular  moment  transitions  change  I 

i 

j 

both  J and  K and  the  lowest  frequency  one  is  at  approximately  A or  | 

5 ^ 

2x10  MHz.  This  is  well  above  the  frequency  range  of  these  experiments.  j 

i 

The  slight  asymmetry  of  the  molecule  splits  each  K>0  level  into  a j 

pair  of  levels.  The  average  frequencies  of  the  A type  transitions  be- 
tween these  split  levels,  within  a good  approximation,  are  given  by  j 

! 

the  average  rotational  constant  (B+C)/2,  i.e.  (B+C) , 2(B+C),  3(B+C), 

etc.  More  accurate  average  rotational  constants  can  be  obtained  by  j 

2 5 i 

applying  small  corrections.  Transitions  between  the  two  split  levels 


J 


of  a single  J state  are  also  allowed.  The  asymmetry  splittings  {and 
also  these  transition  frequencies)  are  strongly  J and  K dependent.  For 


1 


r~ 

\ 

I 
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j ' K=1  and  2,  these  splittings  are  given  (B-C)J(J+1)  and  (B-C)^ 

I 

1 (J-1) J ( J+1) ( J+2) / [16 (2A-B-C) ] respectively.  Water  dimer  should  ex- 

hibit, therefore,  a series  of  radio  frequency  transitions  for  K=1 
asymmetry  doublets  in  the  10  to  1000  MHz  range,  and  in  the  0.1 
to  100  MHz  range  for  K=2  asymmetry  doublets.  For  K>2,  the  transi- 
1 tions  between  asymmetry  splittings  will  be  too  small  to  observe  for 

any  reasonable  J value. 

Several  conclusions  about  water  dimer  radio  frequency  and  micro- 
wave  transitions  should  be  noted.  The  strong  J and  K dependence  of 
j the  radio  frequency  transitions  make  quantum  number  assignment  a 

, trivial  process.  None  of  the  A type  transitions  have  a useful  depen- 

, dence  on  the  A rotational  constants  and  structure  determination  must 

^ be  made  without  this  information.  Stark  effect  measurements  on  A 

5 type  transitions  will  depend  dominantly  on  the  A component  of  the 

permanent  dipole  moment  which  is  nearly  parallel  to  the  oxygen-oxygen 
i axis.  Finally,  each  of  the  two  tunneling  states,  and  E , will  ex- 
hibit slightly  different  effective  rotational  constants.  This  can 

2 2 

be  thought  of  as  a vibration-rotation  interaction. 


Transition  Frequency  Data 

The  first  data  consisted  of  a set  of  radio  frequency  transitions 

in  the  range  of  1-10  MHz  which  fit  to  the  fourth  power  of  the  angular 

momentum;  v = c(J-l)  (J)  (J+1)  (J+2) . This  portion  of  the  spectrum 

clearly  arises  from  asymmetry  splittings  for  K=2  levels.  The  coeffi- 

2 

cient  c in  the  above  expression  is  equal  to  [ (B-C)  /16  ( 2A-B-C)  ] +D 

2 7 

where  is  a centrifugal  distortion  constant. 


While  this  combina- 
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tion  of  constants  is  not  particularly  useful  for  structural  informa- 
tion, the  absolute  certainty  of  the  J and  K quantum  number  assignment 
is  extremely  important  for  identifying  pure  rotational  transitions. 

A second  set  of  transitions  between  K=2  asymmetry  doublets  having  a 
substantially  different  asymmetry  parameter  was  also  observed.  As 
previously  described,  these  two  groups  of  transitions  belong  to  two 
separate  tunneling  levels  having  different  effective  rotational  con- 
stants. Ic  is  not  possible  to  assign  which  set  belongs  to  and  E“ 
representations  and  the  transitions  are  therefore  labeled  E^^  and  E2. 

The  data  and  corresponding  values  for  the  c coefficient  for  H^O  dimer 
are  given  in  table  1.  Similar  results  have  been  obtained  for  (0,0)2 
and  this  information  is  included  in  the  table.  The  coefficient  c is 
not  constant  with  J but  shows  some  additional  centrifugal  distortion 
as  indicated  in  figure  3.  This  J dependence  is  a smooth  function  and 
indicates  t.he  expected  smaller  distortion  for  (020)2.  The  magnitude 
of  the  asymmetry  parameter  would,  however,  increase  on  deuteraticn 
for  a rigid  molecule  contrary  to  the  water  dimer  observations.  It  is 
difficult  to  know  how  to  divide  the  origin  of  this  fact,  as  well  as 
the  large  difference  between  c for  Ej^  and  E2  states,  between  the  rota- 
tional constant  and  0^  portions  of  c.  For  the  usual  rigid  asymmetric 
top  molecule,  the  contribution  to  this  coefficient  would  be  very 
small.  However,  the  expected  very  small  value  of  B-C  and  the  large  A 
value  minimize  the  rotational  constant  term  while  the  low  force  constant 
vibrational  modes  of  water  dimer  enhance  D^.  Anomalous  behavior  of 
this  coefficient  has  been  previously  observed  in  HSSH.  Without  K=1 


information  .it  is  not  possible  to  separate  the  two  contributions. 
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Structural  data  is  best  obtained  from  the  pure  rotational  micro- 
wave  transitions.  To  select  the  desired  rotational  transitions  associ- 
ated with  the  E tunnelinq  levels  from  the  rotation-tunneling  transi- 
tions, double  resonance  techniques  were  employed.  The  J=2,  K=2  to 
J=3,  K=2  pure  rotational  transitions  will  be  at  6(B+C)/2,  or  in  the 
36  GHz  range.  The  microwave  transitions  in  this  frequency  region  were  i 

checked  for  double  resonance  connections  with  radio  frequency  transi- 
tions for  J=^2 , K=2  and  J=3,  K=2  levels.  Since  some  of  the  radio  fre-  ^ 

quency  asymmetry  doublets  for  the  J=2,  K=2  and  J=3,  K=2  states  are  in- 
accessible, the  double  resonance  experiments  were  done  in  the  presence 

of  an  electric  field,  E.  The  radio  frequency  transitions  were  observed 

29 

in  the  strong  field  limit  where  the  Stark  effect  is  linear  and  given 
by  -Uj^EMK/ [J  ( J+1 ) ] . The  dipole  moments,  of  the  two  E states  are 

slightly  different  and  sufficiently  large  electric  fields  were  used  to 
separate  the  two  pure  Stark  transitions.  All  of  the  required  double 
resonance  connections  were  found  and  the  two  microwave  pure  rotational 
transitions  were  thus  identified.  All  that  remained  was  to  identify 
which  microvvave  transition  correlated  to  which  set  of  zero  field  asym- 
metry doublet  transitions.  This  information  was  lost  in  the  double 
resonance  experiment  because  the  high  field  pure  Stark  transitions  do 
not  depend  on  the  zero  field  splitting.  However,  the  correlation  was 
easily  established  by  the  fact  that  the  zero  field  microwave  transitions 
have  observalile  K doubling.  By  tracking  the  high  field  microwave  transi- 
tion to  zero  field  and  then  measuring  the  asymmetry  splitting,  one 
obtains  the  correct  xdentity  of  asymmetry  parameter,  microwave  frt?c]uency, 
and  flipole  moment. 

d 
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The  appropriate  Stark  coefficients  also  confirm  the  quantum 
number  assignment  and  pure  rotational  nature  of  the  microwave  transi- 
tions. In  addition,  the  J=3,  K=2  to  J=4,  K=2  transitions  were  ob- 
served at  four  thirds  the  J=2,  K=2  to  J=3,  K=2  frequencies.  Similar 
measurements  have  also  been  made  on  D^O  dimer.  Table  2 lists  average 
frequencies  for  the  K=2  K doublet  microwave  transitions.  Also  listed 
are  corrected  frequencies  which  account  for  the  small  asymmetry  of 
the  K doublets. 

The  J=3  to  J=4  transitions  are  not  precisely  at  four  thirds  the 
J=2  to  J=3  transition  frequencies  and  this  allows  the  evaluation  of  a 
Dj  centrifugal  distortion  coefficient  for  the  and  states  of  water 
dimer.  This  constant  can  then  be  used  to  estimate  the  monomer-monomer 
stretching  frequency  by  approximating  water  dimer  as  a diatomic  mole- 
cule consisting  of  two  monomers.^®  In  this  case,  D = 4 [ (B+C) /2 ] 

u 

j where  w is  the  harmonic  stretching  frequency  and  (B+C)/2  is  used  as 

the  average  rotational  constant.  While  this  model  is  very  simplistic 
it  is  expected  to  be  relatively  accurate.  The  main  deficiency  is  ig- 
noring angular  vibrational  coordinates  with  small  force  constants,  but 
(B+O/2  is  most  sensitive  to  the  separation  of  the  two  oxygen  atoms  and, 
therefore,  the  stretching  coordinate.  Table  3 lists  and  this  vi- 

u 

brational  frequency  for  the  water  dimer  species  studied. 

While  the  zero  field  splittings  of  the  K=3  asymmetry  doublets  are 
too  small  for  direct  observation,  pure  Stark  transitions  can  be  observed 
and  identified  with  K=3  states.  Using  these  transitions  of  known  quantum 
numbers  in  double  resonance  experiments,  the  J=3,  K=3  to  J=4,  K=3  pure 
rotational  transitions  were  assigned.  This  provides  data  to  obtain  the 
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JK 


centrifugal  distortion  coefficient.  0^.^.  is  rather  large,  as  ex 

u K 


pected  for  a molecule  with  low  energy  vibrational  modes  and  a large 
A value.  K=3  states  have  greater  than  50  cm  ^ of  rotational  energy 
about  the  oxygen-oxygen  axis.  is  an  important  parameter  bectiuse 

it  allows  the  prediction  of  low  K rotational  transitions.  This  mini- 
mizes the  experimental  search,  but  more  importantly,  the  accuracy  of 
the  predictions  is  a measure  of  the  validity  of  using  the  rigid  rotor 
plus  centrifugal  distortion  model  to  interpret  water  dimer  spectra. 

The  K=3  pure  rotational  transition  data  are  given  in  table  2.  The  K=2 

and  3 microwave  data  for  (H2^^0)  ^ {D2O)  ^ were  used  to  predict 

1 8 

similar  transition  frequencies  for  (H2  0)2-  Using  the  limited  supply 

1 8 

of  0 enriched  water  these  transitions  were  observed.  Because  of 
sample  limitations  only  low  resolution  noise  broadened  spectra  were  ob- 
tained and  double  resonance  verification  of  assignments  was  not 
feasible.  The  observed  frequencies  are  listed  in  table  2.  Assictnment 
uncertainties  exist  for  the  K=2,  E_  transitions  and  the  K=3,  E,  and  E-, 
designation  could  be  reversed  producing  only  changes  in  These  un- 

certainties do  not  affect  use  of  these  data  later  in  this  paper. 

Only  or.e  possible  ambiguity  remains  in  the  assignment  of  K=2  and 
3 microwave  data.  For  K=2  transitions  it  was  proven  that  the  lower 
frequency  microwave  data  belongs  to  E^  states  while  the  higher  frequency 
transitions  are  in  E2  states.  The  assignment  of  K=3  data  in  table  two 
assumes  the  same  order.  Reversing  the  order  does  not  yield  D values 
of  particularly  unreasonable  magnitude  but  these  constants  do  predict 
unreasonable  rotational  constants,  (B+C)qq/2,  for  the  J=0,  K=0  rota- 
tionless state.  Using  the  assignment  in  table  2 (B+C)pp/2  for  the 


and  E2  states  are  very  similar,  differing  by  about  2.5  MHz.  This  is 

the  expected  behavior  of  non-rigid  molecules  and  is  the  observed  case 
2 2 3 

xn  ammonia  and  hydrogen  fluoride  dimer.  The  alternate  choice 
yields  constants  over  40  MHz  apart.  This  latter  result  is  not  accept- 
able. Absolute  confirmation  of  E state  assignments  can  be  obtained  by 
the  observation  of  microwave  transitions  for  a third  K state.  For 

several  reasons,  the  most  desirable  data  is  for  K=0  states. 

19 

The  second  order  Stark  effect  of  the  K=0  transitions  requires 
greater  focusing  field  voltages  which  significantly  reduce  the  signal 
to  noise  ratio  of  the  spectrometer.  However,  the  prediction  from  the 
higher  J,  K transitions  proved  to  be  accurate  to  within  one  half 
megahertz  for  the  J=0,  K=0  to  J=l,  K=0  transition,  for  the  E^  state. 
There  are  several  rotation-tunneling  transitions  in  this  frequency 
region  and  careful  confirmation  of  this  assignment  was  necessary.  The 
Stark  effect  of  this  transition  was  accurately  measured  and  confirms 
the  assignment.  This  is  a particularly  rigorous  test  since  seco.id 
order  Stark  effect  requires  that  both  the  dipole  moment  matrix  ele- 
ments and  the  transition  energies  be  correct.  The  J=l,  K=0  to  J=2, 

K=0  transition  was  also  observed  with  a centrifugal  distortion  para- 
meter, Dj,  identical  to  the  higher  J,  K determination.  The  Stark  effect 
of  this  transition  was  also  commensurate  with  the  assignment.  Finally, 
the  J=l,  K-0,  M=0  to  J=l,  K=0,  M=il  pure  Stark  transition  was  observed 
with  appropriate  dipole  moment  and  energy  denominator.  This  transition 
was  also  observed  to  have  a double  resonance  connection  with  the  J=0 , 

K=0  to  J=l,  K=0  microwave  transition. 
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The  K=0  data  are  also  listed  in  table  2.  Only  transitions 
were  observed.  The  transitions  are  masked  by  a stronger  transi- 
tions arising  from  rotation-tunneling  transitions.  All  microwave 
transitions  in  the  12  and  24  GHz  regions  which  have  been  attributed 
to  rotational-tunneling  transitions  have  exhibited  incorrect  Stark 
effect  fov  pure  rotational  transitions.  The  accidental  overlap  of  a 
rotation- tunneling  and  a pure  rotation  transition  is  indicative  of 
the  complexity  of  the  water  dimer  spectrum.  At  the  instrument  sen- 
sitivity required  for  the  K=0  observations,  rotation-tunneling 
transition.s  are  extremely  dense  and  the  overlap  in  question  is  not 
an  improbable  occurrence.  It  has  been  extremely  fortunate  that 
molecular  beam  electric  resonance  spectroscopy  can  be  used  to  easily 
obsein/-e  radio  frequency  transitions,  including  pure  Stark  transitions, 
and  rake  [.recise  dipole  moment  measurements.  This  ability,  along 
with  straightforward  double  resonance  capabilities  between  radic  and 
ricrowave  frequency  transitions,  has  been  essential  in  sorting  out 
water  dimer  spectra. 

Careful  searches  have  been  made  in  the  24  GHz  region  for  the  J=l, 
K=1  to  J=2,  K=1  transitions  without  success.  This  failure  is,  un- 
fortunately, to  be  expected  due  to  the  presence  of  numerous  lines  in 
this  frequency  range  which  are  not  pure  rotational  in  nature  and  the 
large  uncertainty  in  the  asymmetry  doubling.  The  asymmetry  doubling, 
which  could  range  from  0 to  200  MHz,  not  only  makes  the  search  problerr, 
difficult  but  also  determines  the  Stark  effect  and,  therefore,  de- 
termines m.  iny  necessary  experimental  parameters.  Much  effort  was  ex- 


pendcc  in  the  K=1  searches  since  the  asymmetry  splitting  would  provide 
a value  for  B-C  and  very  useful  structural  information. 
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The  K=2  and  3 microwave  transition  frequencies  in  tables  2 have 
been  fit  to  average  rotational  constants,  [(B+C)/2)qq,  for  the  rota- 
tionless J=0,  K=0  state,  D^,  and  D using  the  expression: 

J J K 

V (J,K-*J+1  ,K)  = [ (B+C)-2K^Djj^]  (J+1) -4Dj  (J+1) 
for  each  tunneling  state  of  each  isotopic  species.  The  results,  along 
with  calculated  vibrational  frequencies  are  given  in  table  3.  The  most 
complete  data  set  exists  for  the  state  of  2-  Here  it  is 

possible  to  do  a least  squares  fit  to  1(B+C)/2]qq,  D^,  and  The 

results  are  given  in  table  4.  The  excellent  agreement  between  table 
3 and  4 is  an  additional  confirmation  of  assignments  and  data  treatment. 
Since  the  data  is  the  most  complete  and  reliable,  it  will  be  used 
in  the  following  structure  calculations. 

Stark  Effect  Data 

The  Stark  effect  has  been  discussed  in  the  previous  section  as  an 
aid  to  spectroscopic  assignment.  Here  it  will  be  presented  in  terms 
of  electric  dipole  moment  measurement.  The  electric  dipole  moment  is 
a molecular  property  of  intrinsic  interest,  but  in  water  dimer  the 
dipole  moment  also  plays  an  important  role  in  the  structure  determina- 
tion. This  situation  occurs  because  the  dominant  contribution  to  the 
water  dimer  moment  is  the  vector  sum  of  the  two  monomer  moinents.  There- 
fore, the  relative  orientation  of  the  monomer  units  is  strongly  re- 
flected by  the  dimer  electric  dipole  moment.  In  the  next  section  the 
dipole  moment  data  to  be  presented  here  will  be  combined  with  the  ro- 
tational constant  information  to  obtain  the  structural  parameters  for 
water  dinei . 
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As  mentioned  previously,  pure  Stark  transitions  have  been 

29 

studied  for  states  with  K^2.  Linear  Stark  effect  occurs  when  the 
Stark  energy  is  much  larger  than  the  zero  field  asymmetry  splitting. 

For  K^2  this  condition  is  readily  achieved.  The  Stark  energy  is 
given  by  -u  EMK/ [ J (J+1) ] , where  p is  the  dipole  moment  component 
along  the  A axis.  The  transitions  in  question  are  AM=1  and,  there- 
fore, the  transition  frequency  is  given  by  p EK/[J(J+1)].  For 

this  type  of  transition  has  been  measured  for  K=2,3,4  in 
both  E^  and  E2  states.  Since  the  dipole  moment  is  different  for 
the  two  E states,  two  transitions  will  be  observed  for  sufficiently 
large  fields.  This  will  occur  when  Ap  EK/[J(J+1)]  is  greater  than 

d 

the  linewidth,  were  Ap  is  dipole  moment  difference  between  the  E^^ 
and  E2  states.  Figure  2 shows  this  situation  for  J=2,K=2.  The 
essentially  identical  intensity  of  these  two  transitions  requires  very 
similar  absolute  energy  for  the  E^  and  E2  states.  This  fact,  along 
with  the  observation  of  two  and  only  two  sets  of  transitions,  strongly 
support  the  tunneling  model. 

Accurate  data  was  obtained  for  K=2,3  and  4 states  and  double 
resonance  experiments  were  used  for  the  E^^  and  E2  assignment  as  pre- 
viously described.  The  results,  given  in  table  5,  indicate  significant 
dependence  on  both  K and  E state.  While  this  data  is  unusual,  it  is 
not  surprising  as  it  will  be  shown  that  angle  changes  between  monomer 
units  of  iO.5  degree  will  change  the  moment  by  =.03  Debye.  This 
sensitivity  of  moment  with  angle  change,  along  with  the  35  cm  ^ energy 
difference  between  i'-2  and  K=3  states,  makes  large  changes  in  observed 


moments  the  expected  behavior.  These  measurements  further  show  that 


V 
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the  K dependence  of  the  dipole  moments  is  not  simple  or  monotonic. 

To  observe  any  dependence  on  the  J quantum  number,  J=2 , K=2  and  J=3, 

K=2  moments  were  measured.  The  small  effect  is  shown  in  table  5. 

K=2  and  K=3  dipole  moments  have  also  been  obtained  for  (020)2, 
and  are  listed  in  table  5.  The  moments  for  the  two  tunneling  states 
are  much  closer  together  and  the  K dependence  is  also  much  smaller. 

This  is  consistent  with  (020)2  having  both  smaller  zero  point  vibra- 
tional amplitudes  and  less  centrifugal  distortion.  Because  of  the 
small  difference  between  the  two  (D2O) 2 moments  it  was  not  possible 
to  unambiguously  assign  the  observed  moments  to  tunneling  states. 

The  assignment  in  table  5 is  by  analogy  with  (H^O)  This  assignment 
gives  K and  tunneling  state  dependence  of  the  moment  in  (020>2  which 
is  qualitatively  the  same  as  in  (H2^^0)2* 

The  second  order  Stark  effect  for  K=0  states  is  significantly 

29 

different  from  the  higher  K first  order  effect.  In  principle,  the 
second  order  Stark  coefficients  depend  on  all  three  components  of  the 
dipole  moment  and  also  on  the  A rotational  constant.  However,  in  water 
dimer  u,  and  (B+C)/2  dominate  the  second  order  Stark  effect.  In  an 

a 

attempt  to  obtain  additional  information,  the  Stark  effect  of  three 
second  order  transitions  were  measured  using  long  periods  of  computer 
averaging  to  improve  precision.  The  results  are  given  in  table  6 
for  the  three  transitions  involving  K=0 . The  dipole  moment  values  in 
table  6 were  obtained  from  the  experimental  Stark  coefficients  using 
simple  linear  molecule  second  and  fourth  order  Stark  expressions.^^ 

For  K=0 , this  analysis  is  equivalent  to  a symmetric  top  treatment.  The 
data  is  too  correlated  to  attempt  a full  asymmetric  top  analysis  to 


L 
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cbtain  ^ values.  The  most  useful  approach  is  to  consider 

only  transitions  which  are  allowed  in  the  synmietric  top  limit  .uid  to 

include  only  the  A rotational  constant  in  energy  denominators  for  AK 

transitions . These  assumptions  are  fully  justified  by  tne  precision 

of  the  data.  Using  this  simplified  approach,  the  and  dipole 

moment  components  have  identical  contributions  and  we  can  define 

p = /p,  ^ tp  The  Stark  coefficients  for  the  three  K=0  transitions 
^ DC 

are  given  below. 

J = 0,K=0,M=0-^J=1,K=0,M=0  Au  = [8p^/15  (B+C) +2p^/15A]E  ' 

Si  X 

J=1,K=0,M=0-*J=1,K=0,M=±1  Aw  = [ 3p2/10  (B+C) +p2/5A] 

Si  X 

J = 1 ,K=0  ,M=0-^J=2  ,K=0  ,M=0  Aw  = - [16p2/105  (B+C) +4p2/105A]E^ 

d X 

Note  that  che  two  microwave  transitions,  the  first  and  third  exfiressions , 
have  identical  moment  dependence  and  should  produce  identical  moments 
as  listed  in  table  6.  The  observed  difference  stems  in  part  frc^m  ex- 
perimental uncertainty  and  in  part  from  centrifugal  distortion  effects 

cr.  the  moments.  Since  the  two  data  involving  J=0  and  J=1  energy 

32 

levels  are  more  precise  and  have  less  centrifugal  distortion  t.hese 
can  be  used,  along  with  a calculated  A value,  to  obtain  an  estimate  for 
p . T.his  produces  a value  of  p = .31.15D.  Because  of  uncertain 

1 X 

centrifugal  distortion  effeots this  number  is  not  reliable  but  this  data 
does  allow  reliable  upper  limit  of  0.6D  to  be  placed  on  p^^. 

In  the  following  section  a more  detailed  discussion  of  the  water 
dimier  electric  dipole  moment  will  be  given.  This  will  include  calcula- 
tion of  electrostatic  interactions  between  monomers  and  quantitative 
descriptior.  of  pho  dimer  moment  in  terms  of  the  internal  degrees  of 
freedorr.  of  the  dimer  structure. 
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ETRL’CTURE  CALCULATION 


The  basic  structure  of  water  dimer  can  be  determined  from  experi- 
mental rotational  constants  from  different  isotopic  species  and  from 
electric  dipole  moments.  We  employ  a simple  rigid  rotor  model  to 
calculate  the  dimer  moments-of-inertia  in  the  principal  axis  system 
and  a vector  model,  including  corrections  for  induced  moment  and  charge 
transfer  effects,  to  account  for  the  dimer  electric  dipole  moment. 

The  principal  moments  of  inertia  for  a given  isotopic  species 
of  water  dimer  are  calculated  with  the  assumption  that  the  monomer  sub- 
units have  the  same  geometry  as  the  equilibrium  structure  of  the  free 
water  monomer given  in  Table  7.  Thus  there  are  six  unconstrained 
coordinates  describing  the  orientation  of  the  monomers  in  water  dimer. 
These  are  conveniently  described  by  the  oxygen-oxygen  internuclear 


distance,  Pqq  » and  the  Euler  angles  0^,  ^nd  6^,  Xg(‘*>3  = 0)-  The 


subscripts  refer  to  the  proton  donating  and  proton  accepting  mororaers. 

Ihese  six  coordinates  give  the  orientation  of  the  donor  and  acceptor 

ronorer  molecules  with  respect  to  the  reference  configuration  shown  in 

Fig.  ia.  The  reference  configuration  is  chosen  so  the  axis  of  each 

ronomer  (the  axis  is  assumed  to  point  toward  the  oxygen  atom)  lies 

in  the  positive  Z direction  for  4),  =0j  = 4>  = 6 = 0.  The  and  x 

‘ ^d  d ^a  a '^d  '^a 

specify  the  rotation  of  the  monomers  about  the  axis  with  x^^  = 0 
placing  th"'  donor  protons  in  the  XZ  plane  and  x =0  placing  the  acceptor 
protons  in  the  YZ  plane.  One  of  the  azmuthal  angles  can  be  arbi- 
trarily  set  equal  to  y.fto,  fixing  the  C2  axis  of  the  acceptor  to  lie 
in  the  X/  plane.  The  "trans-1 inear"  structure  predicted  by  calcula- 
tions is  shown  in  Fig.  4b.  This  structure  has  an  XZ  plane  of  symmetry. 
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requirinq  X-,  ~ ~ 0 ^nd  a linear  hydrogen  bond,  fixing  0^  = 

- y a,  where  a is  the  monomer  bond  angle. 

To  obtain  initial  structural  parameters,  the  " trans-1 inear " struc- 
ture of  Fig.  4b  was  used.  This  initial  assumption,  which  will  be 

shown  to  be  valid  below,  can  be  based  on  evidence  from  the  structure 
2 4 

of  ice,  ^ water  dimer  calculations,  and  experimental  and  theoretical 

3 34 

results  on  other  small  hydrogen  bonded  systems.  ' The  most  important 

constraint:  is  the  linear  hydrogen  bond.  This  configuration  minimizes 

repulsions  between  non-hydrogen  bonded  protons  and  is  a reasonable 
35 

first  choice.  The  experimental  evidence  in  its  favor  will  also  be 

B't'C 

presented  below.  From  the  values  for  two  different  isotopic 

species,  R and  0 can  be  calculated.  The  results  are  given  in 
Table  8.  From  the  average  of  the  1112^^0)2/(02^^0)2  and  (H2^^0)„/ 
(D2^^0)2  results,  = 2.98  A and  0^  = 58.8°.  The  (H2^^0) 2 and  (H2^®0) 

results  give  a much  larger  value  for  0 , 82.2°.  However,  relative 

Si 

B+C  16  18 

differences  in  — j—  for  (H2  0)2  and  (H2  0)  depend  primarily  on  R^^ 

and  the  mass  change,  and  only  slightly  on  the  orientations  of  the  mono- 

B+C 

mers.  An  f.rror  of  10  MHz  in  , caused  by  model  errors  (since  the 

experimer.tal  error  is  alm.ost  certainly  a negligible  factor  in  these 

18  16 

calculations),  could  cause  an  error  of  30°  in  0 if  (H_  0)_/(H_  0) „ 

d 2.  2.  2.  2. 

data  are  used  (see  appendix  I) . The  corresponding  error  in  C from 
(H2O)  2/ (020;  2 data  would  only  be  3°.  Thus  the  best  results  shoulil  be 
obtained  from  (020)2/(1120)2  data,  and  simultaneous  fits  of  2 

and  (H2^^*0)2  data  to  calculate  geometries  will  not  work  well  witli  a 
rigid  rotor  model. 


-25- 


To  remove  the  constraint  of  a linear  hydrogen  bond,  it  is  necessary 

B'f’C 

to  use  a third  experimental  datum  in  addition  to  for  two  isotopic 

B+C 

species.  As  discussed  above,  attempts  to  simultaneously  fit  for 

2 2 likely  to  magnify  model  inaccuracies, 


and,  in  fact,  no  structure  would  fit  all  three  values  of 


B+C 


to  better 


than  about  10  MHz,  Fortunately,  the  A-component  of  the  electric  dipole 
moment,  depends  strongly  on  0^  and  0^  and  can  be  used  to  determine 

0,.  It  is  readily  shown  that  y is  given  by  equation  1: 

C dL 


“a  ■ '"ind’a  * '*"cT'a' 


(1) 
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U is  the  electric  dipole  moment  of  the  free  monomer,  1.855D,"'" 
the  dipole  moment  enhancement  caused  by  transfer  of  electronic  charge 
from  the  proton  acceptor  to  the  proton  donor  molecule.  the 

vector  sum  of  the  dipoles  induced  on  each  monomer  by  the  electric  field 
caused  by  t.ie  charge  distribution  of  the  other  monomer.  By  using  a 
nultipole  expansion  for  the  electric  field,  the  dipole  and  quadrupole-in- 
duced  dipole  m.oments  can  be  calculated^^  and  the  results  are  given  in 
table  9.  The  values  of  [y.  for  the  structure  given  below  is  0.42D. 

Since  the  nultipole  expansion  may  not  rapidly  converge  at  the  water  dimer 
internuclea  i:  separation,  and  because  the  charge-transfer  term,  although 
most  likely  sm,all,  may  be  significant,  we  also  obtained  ^ initio  values 
of  the  total  enhancement  of  the  dimer  dipole  moment  due  to  charge-transfer 


and  induced  moment  effects. 


38 


These  results,  found  with  a minimal  basis 


set  are  giver,  in  Table  9,  along  with  the  electrostatic  calculation  results, 

4 

Other  a_b  ini  tio  calculations  with  much  larger  basis  sets  report  total 
dipole  moment  enhancements  in  the  range  of  0.3-0.4D.  All  of  these  results 
predict  a dipole  moment  enhancement  of  about  0.4  i .1  debye.  We  have. 


I 
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therefore,  chosen  to  use  the  electrostatic  calculation  results,  0.42D, 
with  the  very  conservative  uncertainty  limits  of  +0.3D.  With  this 
result  Eq.  1 becomes 

M = u [cos0  + cos6j]  + 0.42  ± 0.30  (2) 

a ad 

The  moment  enhancement  term  has  some  dependence  on  9^  and  6^  and  Eq . 2 

is  valid  for  9^  s 60°  and  6^  = -50°.  In  further  support  of  Eq.  2,  the 

enhancement  of  u in  (HF)_,  over  that  of  the  vector  sum  of  monomer 
^ 2. 

components  is  0.56D,  Thus  Eq.  2 is  quite  realistic. 

B+C 

Using  Eq.  2 as  a third  constraint,  along  with  — for  two  isotopic 
species,  and  0^  can  be  determined.  The  results  are  presented 

in  Table  10.  and  0^  are  almost  unchanged  from  Table  8,  and  is 

within  2°  of  a linear  hydrogen  bond  (6^  = -52.26°  for  a linear  hydrogen 
bond)  . 

The  remaining  constraints  on  the  dimer  geometry,  x ~ Xa  ~ ~ 

are  more  difficult  to  assess.  Xg  and  X(^  give  the  rotation  of  the  mono- 
mers around  their  axes,  and  for  a linear  hydrogen  bond,  gives  the 
torsional  angle  of  in  Fig.  4b,  out  of  the  plane  of  the  drawing. 

Thus  It  is  reasonable,  and  the  calculations  given  in  Table  11  show’  that 
3-t-C 

— which  depends  strongly  on  Rqq»  9^  and  0^,  depends  weakly  on  Xg’ 
and  * The  good  aspect  of  this  is  that  our  values  of  6 and  0, 

are  almost  independent  of  our  choice  Xg  = = 0 . However, 

since  I is  affected  by  X^,  Xj  and  <p,  only  slightly  (the  enhanceir.ent 
term  of  Eq . 2 depends  on  them),  and  at  least  for  the  isotopic  species 
studied  so  far,  is  only  weakly  dependent  on  them,  it  is  necessary  to 

use  otiicr  c.:per  imonta  1 parameters  to  determine  them. 


I 


A 
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The  information  bearing  on  these  coordinates  is  the  upper  limit  ] 

of  0.6D  on  , which  depends  on  <1)^,  as  well  as  0^  and  0^,  and  the 
K=2  doubling  constants.  Table  1,  which  are  strongly  affected  by  x > 
and  Using  a vector  model,  as  in  Eq . 1,  but  dropping  the  en- 

hancement terms  which  will  be  negligible  for  the  purpose  of  this  cal- 
culation , 

~ 2 2 1/2 
y = y [sin  0^  + sin  0,  + 2 sin0  sin0,  cos0,] 

- d Q d Q Q ^ 

] 

Thus  y^  changes  monotonically  from  about  O.ID  to  3.D  as  (J)^  goes  from  i 

0-180°.  Fcr  0 = 58°,  0,  = -51°  and  y ^:0.6D,  a value  of  4),«:  25°  is  I 

Q Q Q 

obtained . 

2 

The  X=2  doubling  constants  depend  on  the  term  (C-B)  /16(2A-B-C)  ! 

which  is  a strong  function  of  as  shown  in  Table  12. 

In  fact,  if  X,  ~ X.J  ~ = Of  C-B  is  less  than  5 MHz  and  water  dimer  i 

a G Q I 

is  nearly  an  accidental  symmetric  top.  If  these  angles  differ  from  ; 

0°,  C-B,  and  (C-B) ^/16 ( 2A-B-C)  rapidly  increase  (Table  12).  As  noted  j 

earlier,  the  experimental  K=2  doubling  constants  do  not  have  the  cor- 
rect isotopic  dependence  for  a rigid  rotor.  Upon  complete  deuteration, 

2 ■ 
(C-B)  /16(2^-B-C)  should  increase  by  about  a factor  of  5 rather  than  : 

the  slight  decrease  observed  experimentally.  Presumably  dynamic  effects 

involving  centrifugal  distortion  (D^  coefficient  in  the  K=2  constant  c) 

and  vibrational  averaging  are  causing  the  deviation  from  rigid-rotor 

behavior.  Assuming  that  such  dynamic  effects  will  be  smallest  in  (020)^, 

the  K=2  doubling  constants  for  that  isotopic  species  can  be  used  to  set 

an  approximate  upper  limit  to  the  deviation  of  x j f X ^rid  from  zero  j 

degrees  (the  near  symmetric  top  case)  . For  c =1.7  )cHz,  an  upper  ] 

2 ! 
limit  of  30'  for  each  angle  is  obtained  (the  corresponding  limit  for  ■ 
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j 

I 


=0.5  kHz  is  only  15°)  . The  situation  becomes  less  clear  if 
and  X,  3re  simultaneously  varied.  If  Xj  ^rid  x ^re  both  increased 

a Q 3 

(C”B) ^ 

from  zero  remains  very  small,  while  if  one  angle  becomes 

positive  and  the  other  negative,  this  quantity  increases  quite  [ 

rapidly.  It  should  be  noted  that  Xj  = 0 for  a linear  hydrogen  bond  ! 

so  it  is  reasonable  that  x^  and  therefore  x^  should  be  < 30°.  The  | 

strongest  statement  based  purely  on  experimental  data  is  that  the  | 

experimental  K=2  doubling  constants  are  consistent  with  values  of 

XjJ  and  = 0,  but  do  not  give  unambiguous  limits  for  Xg  and  x^^  • 

The  s,;ructural  parameters  in  Table  10  represent  an  average 

structure  in  the  sense  that  our  experimental  measurements  are  made  in  I 

the  ground  vibrational  state.  The  vibrational  averages  of  the  measured  | 

quantities  will  vary  from  those  of  the  equilibrium  structure  and  of  1 

1 

different  isotopic  species.  It  is  difficult  to  assess  this  effect  in  ] 


the  absence  of  accurate  vibrational  potentials.  However,  the  dipole 
mcm.ent  component,  u,,  is  constant  within  3%  for  all  isotopic  species 
and  rotational-tunneling  states  measured.  A change  of  only  3°  in  6^ 
or  0^  would  cause  this  variation.  Although  moments  of  inertia  and 
electric  dipole  moments  have  different  operator  dependence,  it  is 
reasonable  that  the  equilibrium  geometry  of  water  dimer  is  within  the 
error  limcits  given  in  Table  10.  By  comparison,  the  much  less  strongly 
bound  ArllCl  (about  0.5  kcal/mole  as  opposed  to  5 kcal/mole  for  (H.,0)^) 

41  4 

complex  shows  a nearly  25%  increase  in  dipole  moment  upon  deuteration.  ' 
The  Ar-Cl-H  angle  calculated  from  the  dipole  moment  is  47.5°  for  ArHCl, 
33.5°  for  ArDCl  and  the  equilibrium  value  of  this  angle  was  thought  to 


! 

i 


i 


1 

be  t 0°.  i 
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In  summary,  we  find  = 2.976(10)A,  0^  = 58°(6)  and  6^  = -51° (6) 
from  the  average  of  the  two  rows  of  Table  10.  From  the  perpendicular 
dipole  moment  analysis,  the  torsional  angle  < 30°.  The  assumption 
that  the  molecule  has  a plane  of  symmetry  (Fig.  4),  is 

consistent  with  K=2  doubling  constant  data,  although  rigorous  limits 
on  Xg  Xj3  cannot  be  determined  from  the  isotopic  speices  studied  to 

date  (partially  deuterated  species  can  resolve  this  problem)  . If 
is  assumed  to  be  zero,  the  above  value  for  9^  gives  the  hydrogen  bond  as 
linear  witiiin  1°  ± 6°. 

DISCUSSION 

4 

A large  number  of  ab  initio  calculations  of  the  structure  of 
water  dimer  have  been  performed,  with  basis  sets  of  varying  size,  and 
in  a few  cases,  with  configuration  interaction.  In  general,"  the 

calculatiors  find  water  dimer  to  have  the  trans  configuration  with  a 
plane  of  symmetry  and  a linear  hydrogen  bond.  The  exact  values  of  the 
ecui librium  geometrical  coordinates  vary  from  calculation  to  calculation. 

O 

Calculations  with  large  basis  sets  have  generally  given  = 3.00  A, 

0^  = 25-58°  and  a linear  hydrogen  bond,  where  this  was  not  a constraint. 

a 

O 

The  agreemeiit  with  the  experimental  R^^  value  of  2.9  8(1) A is  excellent 
and  the  linear  hydrogen  bond  is  in  good  agreement  with  our  value  of  0^  = 

-51(6)°  since  0^  = -52.26°  ( and  x^  0)  gives  a linear  hydrogen  bond  in 
our  coordinate  system.  The  calculated  values  of  0^  are  only  in  moderate 
agreement  with  the  experimental  result  of  58(6)°,  with  the  calculated 
values  generally  being  smaller  than  this  (note  ref.  4m,  however).  Our 
results  concerning  the  perpendicular  dipole  moment  and  K=2  doubling 

] 

J 

r 
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constants  are  compatible  with  the  trans  structure  with  a symmetry 
plane  (tj)^  " Xg  = = 0)  • We  have  a reasonably  good  limit  on  <|)^, 

^ 30°,  from  The  results  bearing  on  and  are  ambiguous,  but 

if  the  hydrogen  bond  is  assumed  to  be  linear,  x.  =0,  the  K=2  dorb- 
ling  constants  suggest  that  x^  30°.  In  summary,  we  feel  that  the 
results  of  the  more  extensive  ^ initio  calculations  are  in  good 
agreement  with  experiment,  except  perhaps  for  0 . 

Table  13  shows  two  interesting  trends  in  this  study  of  (H^O) ^ 

and  the  eailier  study  on  (HF)^  of  Dyke,  Howard  and  Klemperer.^  In 

both  of  these  molecules,  the  heavy  atom  internuclear  distance  is 

considerably  longer  in  the  dimer  than  in  the  corresponding  solid 
4 4 

phase.  Frank  and  others  have  discussed  the  idea  that  cooperative 
phenomena  ai'e  important  in  hydrogen  bonding.  The  results  in  Table  13 
imply  that  hydrogen  bonds  in  the  solid  phase  are  shorter  and  stronger 
than  in  the  dimer.  Apparently  when  a hydrogen  bond  is  formed,  charge 
transfer  occurs  in  such  a way  as  to  enhance  the  formation  of  further 
hydrogen  bonds.  Presumably  negative  charge  is  transferred  to  the 
electronegative  atom  of  the  proton  donor  and  positive  charge  to  the 
hydrogens  of  the  proton  acceptor,  strengthening  the  next  hydrogen  bond. 

A secor d important  point  in  Table  13  is  the  nearly  tetrahedral 
directivity  of  the  proton  acceptor.  It  is  not  at  all  obvious  that  hy- 
drogen bonds,  which  are  weak  compared  to  a normal  covalent  bond,  should 

45 

show  this  di rect ’ onal i ty  . For  example,  electron  density  raps  do  not 
show  prominei'.t  lone  pairs  in  ll^O  or  in  HF,  which  has  cylindrical 
symmetry  an  its  charge  distribution. 


Nevertheless,  unless  this  result 
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is  accidental,  covalent  bonding  apparently  plays  an  important  role 
in  the  hydrogen  bond,  and  in  particular,  in  its  stereochemistry.  In 
regard  to  tiie  important  question  of  the  change  in  orbital  hybridiza- 
tion on  hydrogen  bonding,  i.e.,  ice  has  local  tetrahedral  structure, 
but  the  H2O  bond  angle  is  only  104.52°,  our  calculations  are  not 
sufficiently  accurate  for  a definite  answer.  Studies  of  the  partially 
deuterated  water  dimer,  along  with  a more  refined  model,  may  prove 
to  be  successful  in  determining  this  question. 
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Appendix 

To  shew  that  (B+O/2  for  ^ and  does  not  give 

reliable  parameters,  consider  expanding  (B+C)/2  in  a Taylor  series 
near  the  geometry  given  in  tables  8 and  10.  Retaining  only  linear 
term.s  gives: 


A (B+C) 

ie/2  = 

-4013AR^^  + 

4.O5A0^  + 
d 

4.18A0 

A (B+C) 

00 

N) 

II 

-3631AR^^  + 

3.3OA0,  + 
d 

3.42A0 

A(B+C)^/2  = -3454AR__  + 6.47A0,  + 6.31A0 

D UU  Ct  Si 

where  the  coefficients  are  calculated  derivatives  with  units  of  MHz/A 


and  MHz/decree.  The  derivatives  for  arid  are  all  less  t.han 

0.05  MHz/degree  and  are  neglected.  The  subscripts  16,  18  and  D refer 

to  (H2^^0)2»  (H2^^0)2  and  (D2^^0) 2 respectively.  The  linear  depemdence 

16  18 

of  the  {H2  O) 2 and  (H2  0) 2 is  better  shown  by  dividing  each  equation 

by  the  coefficient  of 


[A (B+C) 

^g/2]/4013  = 

-^^00 

+ 

1.01 

X 

lO" 

^A0  . 
d 

+ 

1.04 

X 

10" 

3a0 

[A (B+C) 

^g/2]/3631  - 

-"^^00 

+ 

0.91 

X 

10" 

^A0  , 
d 

+ 

0.94 

X 

10" 

^A0 

[A (B+C) 

p/2]/3454  = 

- "^00 

+ 

1.87 

X 

lO" 

^A0  , 
d 

+ 

1.83 

X 

lO” 

^A0 

Thus  the  coefficients  of  0^  and  0^  in  the  first  two  of  these  equations 
differ  by  only  10%,  while  the  equation  for  (020)2  has  coefficients 
differing  b/  a factor  of  two  from  the  (H2O) 2 equations.  It  is  easily 
shown  From  these  equations  that  if  the  (B+C)/2  values  are  in  errc"  by 
10  MHz,  e.g.  because  of  zero  point  vibrational  errors,  the  angles  cal- 
culated fron  the  2 (H2^^0)2  equations  could  be  in  error  by 

30°.  If  (020)2  and  either  of  the  (H20)2  equations  are  used,  an  error 
or  only  3°  would  result  from  a 10  MHz  model  error.  (B+O/2  actually  dc- 


pkinds  on  the  dimer  coordinates  in  a non-linear  fashion,  but  it  is 
clear  that  the  (112^^0)2  and  {H2'*’^0)2  data  are  not  sufficiently  inde- 
pendent to  give  reliable  structural  information  in  a simultaneous 
fit.  A more  accurate  model,  including  vibrational  averaging,  will  be 
necessary  to  use  the  (H2^^0) 2 and  (H2^®0) 2 data  together. 
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Table  1.  Radio  Frequency  Transitions 


J 

v„  (MHz) 
^1 

Cp  (KHz)^ 
^'1 

V,,  (MHz) 
^2 

c ( KH  z ) ^ 
^2 

3^ 

0.26 

2.2 

4*^ 

0.26 

0.74 

0.80 

2.2 

5 

0.51 

0.61 

1 .88 

2.24 

6 

0.98 

0 .58 

3.76 

2.24 

7 

1 . 76 

0 .582 

6 . 72 

2 .22 

8 

2 . 88 

0.571 

11.20 

2.222 

9 

4 . 50 

0.56  8 

17.52 

2.212 

10 

6 .72 

0 . 566 

26.20 

2.206 

11 

9.64 

0.562 

37.71 

2.194 

j’ 


_b 

D 

0 ,42 

0.50 

1.40 

1.67 

6 

0 . 76 

0 .45 

2 . 90 

1.7; 

7 

] . 40 

0.463 

5 . 20 

3 . 72 

8 

2 . 30 

0.456 

8.65 

1 . 72 

9 

3.6  0 

0.455 

1 3.5  9 

1.731 

3 0 

5.43 

0.457 

20.30 

1 . 70‘* 

11 

7 . 82 

0.4  56 

29.29 

1.707 

32 

1 1 .05 

0.460 

a . c = 

b.  Honr 

[ (H-C) ^/16 (2A-B-C) 
indor  errors  greatly 

D^l. 

decrease  the 

accuracy  of  c for 

very  lew 

‘‘r'^cuenc';  transit  ions  . 
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Table  2.  Data  from  Microwave  Transitions  J , K-*- J+1 , K“^ 


Ave.Exp. 

Corrected 

Species 

J,K 

Freq . (MHz ) 

Freq.(MHz ) 

0,0 

12 

321.00 

1,0 

24 

640.88 

2,2 

36 

863.51 

36 

863.46 

3,2 

49 

146.13 

49 

146.07 

3,3 

48 

985.55 

i 

: (h/<>0)2E2 

2,2 

36 

928.57 

36 

928.45  i 

3,2 

49 

232.47 

49 

232.23  j 

3,3 

49 

205.54 

1 

1 

1 

(D2^^0) 

2,2 

32 

549.20 

32 

1 

549.18 

3,2 

43 

395.24 

43 

395.19 

3,3 

43 

325 . 70 

; 

2,2 

32 

564.58 

32 

564.50 

3,2 

43 

415.58 

43 

415.41  j 

3,3 

43 

366.63 

1 

1 

1 ft 

(H^  0) 

2,2 

33 

333.0 

c 

3,2 

44 

440.2 

c 

3,3 

44 

297.2 

2,2 

33 

384.5 

c 

3,2 

44 

507.8 

c 

1 

3,3 

44 

301.5 

Experimental  uncertainties  are  in  least  significant  figure  quoted. 

Corrections  applied  to  the  average  of  the  two  K doublet  frequencies 

generate  frecpaencies  dependent  only  on  (B+O/2,  see  ref.  25.  These 

corrections  are  vanishingly  small  for  all  but  K=2  transitions. 

Insufficient  data  are  available  for  corrections  to  be  made. 

1 8 

There  are  uncertainties  in  the  assignment  of  F'2  O frequencies. 
These  data  arr  included  only  for  completeness  and  assignment  may  be 
in  error. 


Table  3.  Properties  Calculated  from  K=2,3  Microwave  Frequencies 


Species 

[ (D+C)/2)]qq  (Milz) 

D,  (KHz) 

kj 

00  (cm 

Djj^(MHz) 

"1 

6 

160.8 

46.7 

149 

4.01 

z z 

^2 

6 

158.3 

50.1 

143 

0.67 

E, 

5 

432.4 

33.1 

147 

1.74 

Z Z 

^2 

5 

4 32.9 

35.0 

143 

1.22 

'•1" 

5 

570.4 

34.1 

149 

3.57 

a.  D given  here  depends  on  £'  assignment,  see  text.  E„  assignments  are 
j K 18^  ^ 

uncertain  tor  (H2  0)^/  see  text  and  table  2. 


L 


1 
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Table  4.  Least  Squares  fit  Properties  of  2 


[(BtC)/2]oo 

- 

6 160.7(1) 

MHz 

°JK 

= 

4.01(1) 

MHz 

= 

44  (4) 

KHz 

f 
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Table  6 . K=0  Dipole  Moment  Data  for  2 


J , K , M-*- J ' , K , M ' 

Stark  Coefficient^ 

C 

2 

(Debye) 

0 , 0 , 0 -*•  1 , 0 , 0 

[8m  ^/15 (B+C) +2p^/15A]E^ 

7.66301 (20) xl0“^ 

2.64307 (5) 

1,0,0-1,0,±1 

[3u^^/10  (B+C)+uf/5A]E^ 

4.31197(30)xl0"^ 

2.64354  (10) 

1 , 0 , 0--2 , 0 , 0 

-[16m  ^/105(B+C)+4m^/105A]E^ 
a i 

-2.19409 (80) xlO"^ 

2.64588(75) 

a.  See  text  for  definitions  and  assumptions. 

2 4 

b.  Data  was  fit  to  the  form  v = C_  + C_E  + C.E  . C.  was  calculated  and  C 

0 2 4 4 

and  C2  determined.  C2  is  given  in  cgs  units. 

c.  This  moment  obtained  from  C2  and  the  first  term  in  the  listed  "Stark 
Coefficient".  Values  are  relative  to  u(OCS)  = 0.71521  (see  ref.  31a). 
Numbers  in  parentheses  are  one  standard  deviation  in  the  precisicn; 
absolute  accuracy  is  ±0.02%. 


i 

i 
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Table  7.  Water  Monomer  Properties  Used  in  Obtaining 
Water  Dimer  Structure  Parameters 


0-H  equilibrium  bond  length 

r 

e 

= 

0.9572  A 

Re  f 
33 

H-O-H  equilibrium  angle 

“e 

= 

104.52° 

33 

electric  dipole  moment 

= 

1.855  D 

36 

polarizability  component 

“xx 

= 

1.452  A-* 

39 

ayy 

= 

° 3 

1.226  A-^ 

39 

“zz 

= 

o 3 

1.651  A-^ 

39 

electric  cjuadrupole  moment 

4> 

XX 

= 

0.13  X 10  esu-cm^ 

40 

YY 

-2.63  X 10  esu-cm^ 

40 

zz 

= 

— 9 ft  T 

2.50  X 10  esu-cm 

40 
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B+r* 

Table  8.  and  0^  from  (~^)  for  two  isotopic  species. 

B+C 

The  (— ^)  values  are  for  J=0,  K=0  and,  for  the 

sake  of  definiteness,  the  tunneling  state. 

For  reasons  discussed  in  the  text,  the  2 

1 8 

and  (H2  0)2  results  are  considered  to  be  un- 

reliable . 


Isotopic  Species 

^00 

0 

(A) 

8 

a 

(Degrees 

(Hj“0)2 

and 

2 

2.980 

60.7 

1 Q 

(H2  0)2 

and 

(D2^%2 

2.973 

56.8 

(H2^^0) 2 

and 

(H2^®0)2 

3.000 

82.2 

3.000 


82.2 
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Table  9.  Electrostatic  calculation  of  induced  dipole 
moment  at  several  geometries. 


o 

o 

a: 

®d 

0 

a 

^^ind^  a 

^^ind'a  ^ ^^*^CT'a 

(ab  initio . ref . 38i 

2.9  8 

-52.2 

58 

0.42 

0.53 

2.9  8 

-57.2 

63 

0.41 

0.52 

2.98 


-47.2 


53 


0.42 


0.53 
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Table  10.  ^^qO'  ®a'  ~2~ 

error  limits  in  parentheses  are  calculated  by 
assuming  the  dipole  moment  enhancement  to 
have  an  error  of  ±0.30. 
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Tdble  11 . Variation  of  for  (H2^^0)  ^ with  angular  co- 

O 

ordinates.  = 2.980  A.  Comparison  of  each 

B+C 

geometry  with  the  first  shows  that  depends 

strongly  on  0^  and  0^,  but  considerably  less 

soon  if  j , X T 3nd  x • 
d ^d  '^a 


^d 

Degrees 

r. 

d 

Degrees 

Xa 

Degrees 

0 

a 

Degrees 

^a 

Degrees 

B+C 

2 

NUiz 

0 

-52 . 26 

0 

58 

0 

6 

150.9 

0 

-52.26 

0 

0 

0 

6 

020.2 

0 

0 

0 

58 

0 

6 

274.5 

0 

-52.26 

0 

58 

90 

6 

122  . L 

0 

-52.26 

90 

58 

0 

6 

180.  .5 

180 

-52.26 

0 

58 

0 

6 

160.3 

0 

0 

0 

0.0124 

1 

0 

0 

2.33 

i 0 

30 

0 

2.12 

0 

0 

30 

1.65 

0 

30 

30 

0.00771 

0 


30 


-30 


5.20 
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Figure  Captions 


Figure  1: 
Figure  2: 


Figure  3: 
Figure  4 : 


Cross  sectional  views  of  mechanical  microwave  modulator. 

Pure  Stark  transitions  for  the  and  J=2,  K=2  states. 

These  results  were  obtained  using  signal  averaging  with  an 
effective  time  constant  of  100  seconds.  The  linewidth  in- 
cludes unresolved  hyperfine  structure. 

Radio  frequency  K=2  transitions  plotted  as  a function  of  J. 

a.  Reference  configuration  for  water  dimer  structure.  The 
Y axis  points  toward  the  observer  with  the  acceptor  in 
the  YZ  plane  and  the  donor  in  the  XZ  plane. 

b.  Water  dimer  configuration.  0^  and  6^  give  the  ancles 

m.ade  by  the  C2  axis  of  the  acceptor  and  donor  monomers, 

respectively,  with  the  Z (oxygen-oxygen)  axis.  The 

azm.uthal  angle  is  fixed  at  zero  and  then  gi /es 

the  torsional  angle  between  the  monomers,  arid 

specify  the  rotation  about  each  monomer  C_  axis.  ) = 

2 'a 

= 4)d  = 0 specifies  a dimer  with  an  XZ  symmetry  plane. 
A linear  hydrogen  bond  has  and  x^_j  = 0,  wliere  a 

is  the  monomer  bond  angle. 
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